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MAJOR MERGING: THE WAY TO MAKE A MASSIVE, PASSIVE GALAXY 
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ABSTRACT 

We analyze the projected axial ratio distribution, p(b/a), of galaxies that were spectroscopically 
selected from the Sloan Digital Sky Survey (DR6) to have low star formation rates. For these quiescent 
galaxies we find a rather abrupt change in p(b/a) at a stellar mass of ~ 10 M©: at higher masses 
there are hardly any galaxies with b/a < 0.6, implying that essentially none of them have disk-like 
intrinsic shapes and must be spheroidal. This transition mass is ~ 3 — 4 times higher than the 
threshold mass above which quiescent galaxies dominate in number over star-forming galaxies, which 
suggests that these mass scales are unrelated. At masses lower than ~ 10 11 M Q , quiescent galaxies 
show a large range in axial ratios, implying a mix of bulge- and disk-dominated galaxies. Our result 
strongly suggests that major merging is the most important, and perhaps only relevant, evolutionary 
channel to produce massive (> 10 11 Mq), quiescent galaxies, as it inevitably results in spheroids. 
Subject headings: galaxies: elliptical and lenticular, cD — galaxies: formation — galaxies: fundamental 
parameters — galaxies: statistics — galaxies: structure 



1. INTRODUCTION 

Even galaxies with little star formation activity con- 
tinue to evolve, as evidenced by the substantial in- 
crease of their cos mic stellar mass density over the 
past 7 billion yea rs (jBell et al.|[200l IFaber et all 120071 : 
iBrown et al.ll2007h . This must be related to the decreas- 
ing star formation ac tivity over the same period (e.g., 
iLe Floc'h et alJl2005[ ). and the production of such qui- 
esce nt galaxies through the truncation of star formation 
(e.g.. IFaber et al. 2007; Bell et aU200l : the color scatter 
among quiescent galaxies and its e volution are in precis e 
agreement with such a scenario (Ru hland et alJl2009h . 
There are, however, quiescent galaxies at all redshifts 
z < 1.3 that are more massive than the most massive 
star- forming galaxies. This implies that star formation 
in the most massive galaxies was truncated even earlier, 
and/or that mergers play an important role in producing 
massive galaxies. 

Evidence for the early formation of massive galaxies is 
provided by their old stellar populations. However, we 
need to bear in mind that there can be a large differ- 
ence between the age of the stellar population and the 
assembly age, especially if mergers are important, as is 
the case in a hierarchi cal framework for galaxy formation 
(|De Lucia et al.ll2007h . Hence, the number density evolu- 
tion of galaxies is important in constraining their assem- 
bly history. Measuring this is difficult because of its sen- 
sitivity to the luminosity evolution correction, especially 
for massive galaxies at the exponential cut-off of the mass 
function. As a result, there is no c onsensus among the 
currently available measurements (iCimatti et a.1.1 120061 : 
IWake et/al.ll2lMlBrown"et~aIll2007l : ICool et alJl2003 T^ 

Given these difficulties, other observations have been 
used to either directly or indirectly constrain the assem- 
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bly of galaxies. Merging activity among the massive 
galaxy population i s obs e rved (e.g., van Dokkum et al 

2006al |bl: iLin et al 
produce a color- 



1999; Ivan Dokkuml 120051 : iBell et all 
2008), and has been shown to 



magn itude relation that is in agreement with observa- 
tions (jSkelton et al.ll2009f ). However, its cosmological rel- 
evance has always been difficult to determine, given the 
uncertainties in converting observed merger fractions to 
merger rates and the associated growth in mass. An in- 
dependent and indirect indication that massive galaxies 
undergo continuous evolution is provided by the recent 
result that high-redshift quiescent galaxies are substan- 
tially smaller than local g alaxies with the same mass (see, 
Ivan der Wel et all 120081 and references therein) . This 
strongly suggests that m ergers are important (see, e.g., 
Ivan der Wel et ahl 120091 ). and that the assembly of mas- 
sive galaxies is continuing up until the present day. An- 
other indirect, yet powerful, constraint is provided by 
the evolution in the clu stering and halo occupation dis- 
tribution of red galaxies ({White et al.l2007t iConrov et al.l 
120071: IBrown et al.ll2008h : the evolution in the clustering 
strength of red galaxies is slower than expected in the 
absence of merging. 

In this Letter we address the question whether major 
merging is the dominant mechanism for the production 
of very massive, quiescent galaxies. The argument that 
we invoke is simply that major merging generally leads to 
rounder galaxies. An analysis of the shape distribution 
of quiescent galaxies can therefore constrain the impor- 
tance of merging. Since merging among galaxies with 
mass ratios of < 3 is the only known mechanism to pro- 
duce round galaxies (see Section 3 for further discussion) , 
this is a powerful test. The disadvantage of this method, 
compared to those mentioned above, is that no informa- 
tion about the time scale and epoch of galaxy assembly 
c an be inferred. 

IVincent & Ryder] ((20051) and lPadilla fc Straussl (l2008h 
were the first to systematically study the axial ra- 
tio distribution, p(b/a), of a large number of galax- 
ies, selected from the Sloan Digital Sky Survey (SDSS). 
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Fig. 1. — Left: axial ratio distribution, p(b/a), as a function of stellar mass for spectroscopically selected quiescent galaxies from the 
SDSS at 0.04 < z < 0.08. The gray scale represents, normalized to the total number of galaxies in narrow bins of stellar mass, the fraction 
of galaxies with axial ratio b/a. The upper boundaries below which, as a function of mass, 10%, 30%, 50%, 70%, and 90% of galaxies are 
located, are delineated by the red lines. Right: fraction of galaxies with axial ratios b/a smaller than 0.8, 0.6, 0.4, and 0.2, as a function 
of stellar mass. These figures clearly show that at M* > 10 11 Mq, the fraction of galaxies with small b/a decreases rapidly with mass. 
At lower masses, p(b/a) is approximately uniform in the range 0.3 < b/a < 0.9, implying a significant contribution of disks. At higher 
masses, axial ratios are approximately evenly distributed in the range 0.6 < b/a < 0.9, which shows that disks must be rare, and galaxies 
intrinsically round. 



Through a detailed analysis, they infer the intrinsic 
shape distribution and the effect of extinction. Both 
divide the sample into 'elliptical' and 'spiral' galaxies, 
and confirmed that luminous 'elliptical' galaxies are, 
on average, rounder and tri-axial, compared to low- 
luminos ity 'ellipticals', which are more elon gated and 
oblate (jDavies et al.lll985t iFranx et alj 119911). a nd dis- 
play disky isophotes dJorgensen fc Franxl 119941) . This 
phenomenon is not recent: iHolden et alT(|2009( ) showed 
that this trend persists at least out to z ~ 1. Here we 
present a complementary, modified analysis, focusing on 
p(b/a) as a function of stellar mass for quiescent, i.e., 
non-star- forming, galaxies. Because mass-to- light ratios 
are well constrained by broad-band colors for quiescent 
galaxies, stellar mass estimates are robust. This is es- 
sential for our purposes, as we are interested in the most 
massive objects, i.e., those that populate the exponen- 
tial tail of the mass function. Furthermore, as opposed 
to previous studies, we pre-select galaxies independent 
of their photometric properties. Our shape-independent, 
spectroscopic selection criteria circumvent the biases that 
are potentially introduced by selecting galaxies by their 
'morphological' properties, or some pre-defined surface 
brightness profile. 

With this sample, for which we have determined ax- 
ial ratios from our own fits to two-dimensional light dis- 
tributions, we address the following specific questions. 
Are high-mass, quiescent galaxies rounder than low- 
mass quiescent galaxies? If so, is there a mass limit at 
which p(b/a) distinctly changes, and above which disk- 
dominated are completely absent? Such evidence would 
imply that the only evolutionary path to such masses is 



a disk-destroying mechanism, i.e., major merging. 
2. THE SAMPLE 

We select a sample of 17, 480 quiescent galaxies from 
Data Release 6 of the SDSS (| Adelman-McCarthvl[200l ) . 
Our sample includes galaxies at redshifts 0.04 < z < 0.08 
without detectable [Oil] and Ha emission lines. The se- 
lection criter i a are described and motivated in full by 
IGraves et al.l (|2009D ; but as opposed to that work, we do 
not exclude galaxies with a low concentration index and 
gala xies that are fit better by an exponential profile than 
by a ide Vaucouleursl (| 19481 ) profile, because this may ex- 
clude quiescent, yet disk-like galaxies, which are obvi- 
ously relevant for quantifying p(b/ a) of quiescent galax- 
ies. As a consequence, our sample may include galax- 
ies with star formation in an extended disk outside the 
SDSS spectroscopic fiber. This effect, however, does not 
affect our main conclusion that quiescent massive galax- 
ies with prominent disks are extremely rare (see Section 
[3]). Rather, such a bias works in the opposite direction 
in the sense that it would lead to the mistaken inclusion 
of galaxies with large disks. 

The exclusion of all galaxies with emission lines also 
excludes quiescent galaxies with active galactic nuclei. 
Their number, however, is sma ll, and make up a small 
fraction of the population fe.g.. lPasauali et "aill2009f ) that 
is negligible for our purposes. 

The axial rati o s we re obtained as de scribed by 
van der WeleTaT] (pOOST ). Briefly, GALFIT (|Peng et all 



2002f ) is used to determine from the r-band the radii, ax- 



ial rat ios, position ang l es, an d total magnitudes, assum- 
ing a Ide Vaucouleursl 1)19481 ) surface brightness profile. 
We have verified that adopting surface brightness mod- 
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els with a free Sersic index does not lead to a significantly 
different p(b/a). 

The stellar masses are derived with t he simple conver - 
sion from color to mass-to-light ratio (iBell et al.l [20031. 
but are normalized to correspond to the iKroupal (|2001l ) 
stellar initial mass function. The assumed cosmology is 
(n M , n A , h) = (0.3, 0.7, 0.7). 

The sample is complete over the entire redshift range 
0.04 < z < 0.08 down to M, ~4x fO lo M , set by the 
spectroscopic magnitude limit of the SDSS (r — 17.7). 
The SDSS may be incompl ete for low-luminosi ty, low- 
surface brightness galaxies (Blanton et al.ll2"005[ ). which 
could, in addition, de pend on their orientation (see, e.g., 
IQdewahn et al.l Il9 971 ). However, since we are concerned 
with the high-mass end of the galaxy population, this 
does not play a role. Moreover, simulations of images 
with even lower signal-to-ratio than those of the mas- 
sive galaxies analyzed here demonstrate that axial ratio 
measurements from GALFIT are robust and accurate 
(IHolden et al.ll2009lh In summary, the lack of galaxies 
with small b/a, reported in the following section, is not 
in any way compromised by selection effects or measure- 
ment errors. 

3. RESULTS AND DISCUSSION 

In Figure [lja) we show p(b/a) of the 17,480 spectro- 
scopically selected, quiescent galaxies as a function of 
stellar mass. p(b/a) is shown in gray scale, with the per- 
centiles of the cumulative b/a distribution shown as (red) 
lines. FigureUfa) immediately demonstrates that for qui- 
escent galaxies, the projected axial ratio distribution is a 
strong function of stellar mass. In the narrow mass range 
8 x 10 10 < M»/M© < 2 x 10 11 there is a rapid decrease 
in the number of galaxies with small axial ratios. As fur- 
ther illustrated by FigureQlb), above M* - 2 x 10 11 Mq 
quiescent galaxies with b/a < 0.6 are essentially absent. 

This result shows that evolutionary paths that lead to 
quiescent galaxies with stellar mass M* > 2 x 10 11 Mq 
all but exclude the existence, or the survival, of highly 
flattened, disk-like stellar components. As highly flat- 
tened stellar systems are quite common at lower masses, 
in the possible realm of plausible progenitors of high- 
mass galaxies, this result implies the destruction of the 
flattened component in whatever process causes growth 
beyond M* ~2x 10 11 Mq. Therefore, our result that 
essentially all quiescent galaxies with masses larger than 
M* ~2x 10 11 Mq are round strongly suggests that for 
such galaxies major mergers are the dominant, perhaps 
even unique, formation channel. The destruction of a 
stellar disk requires a major merger, i.e., a merger involv- 
ing progenitors with a relatively small mass ratio of at 
most ~ 3, mergers with a larger mass ratio leaving stel- 
lar d isks intact (see, e.g., iBekkil IT998I : iBournaud et al.1 
I2004T ). Moreover, most likely, the progenitors are not 
very gas rich, as this would produce a disky remnant 
(e.g. jNaab et al.ll2006l ). 

It has been suggested that cold flows are responsible 
for th e formation of mas s ive, classical bulges at high red- 
shift toekel et al.l 120091 : iCeverino et all 120091 ) . In this 
scenario, intensely star-for ming 'knots' m erge, forming 
a massive bulge (see also iNoguchil [1999) . However, a 
substantial fraction of the mass (~ 50%) is still pre- 
dicted to reside in a disk. Even at later stages, when 
the gas disk has become stable against fragmentation 



and collapse ( "morphological quenching" ; iMartig et al.1 
2009), the stellar disk remains intact and contains a 
non-negligible fraction of the total mass. In short, al- 
though cold flows plausibly produce quiescent galaxies, 
the end-products will not be uniquely round. Only in the 
case of sufficient merger activity would galaxies become 
spheroidal. 

In passing, we note that the sharp decrease in the 
fraction of very round galaxies (b/a > 0.8) at the very 
highest masses (M > 3 x 10 11 Mq) signifies that such 
high mass galaxies are typically brightest group/cluster 
galaxies, which tend to be slightly more elongated than 
'norm al' massive elliptical galaxies (see, IBernardi et al.l 
120081) . 

As already noted above, at masses lower than M* ~ 
10 11 Mq, quiescent galaxies display a large range in ax- 
ial ratios, which implies that star formation truncation 
mechanisms below 10 11 Mq are often not associated with 
the destruction of the disk. It remains to be tested 
whether p(b/a) of low-mass quiescent galaxies is simi- 
lar to or different from p(b/a) of star-forming galaxies 
in the same mass range. Such an analysis, which is non- 
trivial because of the effects of extinction and color gradi- 
ents, will constrain the degree to which mergers or bulge 
growth regulate star formation at these lower masses. 

Another open question concerns the number and prop- 
erties of massive, star-forming galaxies. Morphological 
studies suggest that a large fraction (20% — 40%) of all 
galaxies more massive than M ~ 10 11 Mq are lat e-type 
galaxies (|van der Well l2008t iBamford et all l2009f ). and 
their high masses of at least some of these objects are con- 
firme d by their rotational velocities (e.g.. lCourteau et al.l 
2007). Yet, the de gree to which such galaxies are disk- 
dominated and should be considered actively star form- 
ing remains to be determined. It would, therefore, be 
premature to conclude that merging is the only way to 
produce a massive galaxy in general, and therefore re- 
strict this proposition to the formation of massive, qui- 
escent galaxies. 

The picture sketched by the axial ratio distribution of 
quiescent galaxies is in agreement with the strong cor- 
relation b etween structure and mass for galaxies in gen- 
eral (e.g. JKauffmann et al.|[20Tjllvan der We]||2008D and 
early-type galaxies in par ticular (e.g.. iCaon et al.l [19931 : 
iGraham fc Guz man 2003): high-mass galaxies are more 
concentrated and have higher Sersic indices than low- 
mass galaxies. These trends are an indirect indication of 
a decreasing importance of disks for galaxies with higher 
masses, although part of this trend is caused by the in- 
crease in Sersic index with galaxy mass among spheroidal 
galaxies. In our sample we see a similar trend: in the 
mass range 10 10 < M*/M < 2 x 10 10 , 41% of the galax- 
ies have Sersic indices n < 3, whereas at higher masses, 
M* > 2 x 10 11 Mq, only 3% have such low Sersic indices. 
We postpone a full exploration of the joint behavior of 
shape and structure as a function of galaxy mass until a 
future paper, but it is encouraging that the apparent ab- 
sence of prominent disks in high-mass, quiescent galaxies 
is reflected in both the Sersic index and p(b/a). 

The authors thank the referee for positive feed- 
back. A.v.d.W. thanks Marijn Franx for helpful dis- 
cussion. H.W.R. thanks Simon White for his insis- 
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tence on eventually 'publishing the last thesis chapter'. 
This Letter makes use of the Sloan Digitial Sky Survey 



(www.sdss.org). 
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